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We explore the dynamics of dark energy models based on a Dirac-Born-Infeld (DBI) tachyonic 
action, studying a range of potentials. We numerically investigate the existence of tracking be- 
haviour and determine the present-day value of the equation of state parameter and its running, 
which are compared with observational bounds. We find that tachyon models have quite similar 
phenomenology to canonical quintessence models. While some potentials can be selected amongst 
many possibilities and fine-tuned to give viable scenarios, there is no apparent advantage in choosing 
a DBI scalar field instead of a Klein-Gordon one. 



PACS numbers: 98.80.Cq 



I. INTRODUCTION 



Most dark energy modelling using scalar fields has 
followed the quintessence paradigm of a slowly rolling 
canonical scalar field. However, there has been increas- 
ing interest in loosening the assumption of a canonical 
kinetic term. In its most general form, this idea is known 
as k-essence 0. A more specific choice is the 'tachyon' 
which can be viewed as a special case of k-essence 
models with Dirac-Born-Infeld (DBI) action Q. This 
kind of scalar field is motivated by string theory as the 
negative-mass mode of the open string perturbative spec- 
trum, though its use in the dark energy sector is primar- 
ily phenomenological. One goal of such studies is to in- 
vestigate whether there are any distinctive signatures of 
non-canonical actions available to be probed by observa- 
tions. For a recent comprehensive review of dark energy 
dynamics, see Ref. 

Tachyon dark energy has been explored by many au- 
thors, for example Refs. 0, B II ll H ES-^ Two pa- 
pers are particularly closely related to the present work. 
Bagla et al. focussed on two specific choices of tachyon 
potential, and carried out numerical analysis of the cos- 
mological evolution in order to constrain them against 
supernova data and the growth rate of large-scale struc- 
ture. Copeland et al. y\ studied a wider range of po- 
tentials, concentrating mainly on analytical inspection of 
attractor behaviour and the critical point structure with- 
out making comparison to specific observations. In this 
paper, we aim to merge some of the positive features of 
each analysis, by studying a wide range of potentials and 
testing them directly against current observational con- 
straints as given in Ref. [l3l |. 

The mechanism of slow rolling is the key ingredient in 
order to get an accelerating evolution driven by a scalar 
field. Since the DBI action can be expanded to match 
the Klein-Gordon one in this regime, one does not ex- 
pect to find radically different features from the tradi- 



tional quintessence models. However, we point out that, 
as compared to canonical quintessence, tachyon models 
require more fine-tuning to agree with observations. This 
is consistent with the properties of the slow-roll corre- 
spondence between the tachyon and an ordinary scalar 
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II. TACHYON DARK ENERGY MODELS 

A. Equations of motion 

We assume a four-dimensional, spatially-flat 
Friedmann-Robertson-Walker Universe filled by dust 
matter (subscript 'm'), radiation ('r') and a minimally 
coupled homogeneous DBI tachyon T with potential 
V{T) and dimension E~-^. For a perfect fiuid n with 
energy density pn and pressure p„, the barotropic 
index is w„ = Pn/ Pn- Each fluid component satisfies a 
continuity equation 



3ffp„(l-Hw„) = 0, 



(1) 



with w-ai = and Wj- = 1/3. Here, a dot is derivation with 
respect to synchronous time and H = a/a is the Hubble 
parameter defined in terms of the scale factor a{t). In the 
following a subscript will denote quantities evaluated 
today (at to), when a{tQ) — ao — I. Defining the critical 
density today as pco = SHq / (SttG) and assuming that 
gravity obeys an Einstein-Hilbert action, the Fricdmann 
equation reads 



(2) 
(3) 
(4) 





= Pm + ffr + 












U{T) 


Qt 





Vi - r2 ' 



(5) 



where = Pn/ Pco, U 
has been rescaled as t 



^ Here we d o not consider the tachyon either as a dark matter 
candidate Illlll2l or as the infiaton field. 



y/PcO, and the time coordinate 
t/Ho so that H today. The 
tachyon is dimensionless in these units. Note that the Qn 
are the densities normalized to the present value of the 
critical density, and are not the density parameters. The 
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present density parameters have values rim.o ~ 0.24 
and n,^o « 8 X lO'^^ 

The tachyon equation of motion is 



T 



i-r2 



3HT 



IT 



= 0, 



(6) 



where U^t = dU/dT. Since pr = -V\/l-t^, the 
barotropic index for the tachyon is 



1, 



(7) 



which can vary only between —1 and in order for the 
action to be well defined. When the scalar field slowly 
rolls down its potential, \T\ <C 1, it behaves like an effec- 
tive cosmological constant, wt ~ — 1. From now on we 
drop the subscript T' on wt- 

Since the cosmological evolution spans many orders of 
magnitude in synchronous time, for numerical work it is 
convenient to switch to the number of e-foldings N = 
Ina/oo = In a as the evolution parameter, so that 



m,oe 



-3N 



(8) 



Derivatives with respect to N will be denoted by primes, 
so that Q = HQ' and Q = H'^{Q" ~ eQ') for any quan- 
tity Q, where e = -H/H^ = -{H^y/{2H^). With the 



compact notation x 
equation becomes 



u 



Qr, the Friedmann 



Vl - xT'^ 

The equation of motion for the tachyon is 
x{T" -eV) , „ ^, , C/,T 



1 - xT'^ 



where 



X 

2x 



•ixT 



U 



= 0, 



ixUT'^ 
Vl - xT'^ 



From Eq. H1U|I one finds 



rj-i/l 



(1 - XT'- 



3xT' + 



(9) 



(10) 



(11) 



(12) 



Mapping between the number of e-foldings and the red- 
shift z = ap/a — 1, we note that at big bang nucleosynthe- 
sis (BBN) -/Vbbn ~ —20 (z ~ 10^), at matter-radiation 
equality iVcq w —8 (z w 3200), and at recombination 
iV,ec«-7 (z«1100). 

As regards the initial conditions at early times for 
the dynamical equations, we can consider two qualita- 
tive cases. In the first one, the scalar field starts rolling 



^ As input in our numerical code we have chosen fJm.o = 0.25. All 
results are unaffected by small changes in the matter density. 



down very slowly, xT'"^ <C 1. The Friedmann equation 
becomes linear in x and, during radiation domination, 
Q ^ U and x Pi Q. In such a picture of the events, 
the initial condition is a radiation-dominated Universe 
in which a dynamical cosmological constant is negligible. 
Later on, qt increases relative to gr, eventually domi- 
nating at low redshift, while w stays negative and varies, 
perhaps non-monotonically, from around — 1 to < 0. 
The parameters of the model and the shape of the po- 
tential can be adjusted so that the actual value of wq is 
compatible with supernovas data. 

In the second case xT''^ < 1, that is, T' is large enough 
to compensate, but not override, x ^ 1. The tachyon 
starts from a dust regime and ends up again with a suit- 
able negative wq. This is precisely the situation studied 
numerically in Ref. ^] for U cx . 

Below, we shall focus mainly on the first case, which 
encodes all the relevant features of the models. 



B. Tracking and creeping regimes 

Quintessence behaviour typically falls into one of two 
classes, named as tracking and creeping in Refs. [mil 
(see also Refs. [23, 0| which refer to the more general 
late-time behaviour in these scenarios as freezing and 
thawing). The same kind of evolution also occurs in the 
tachyon case. 

Tracking behaviour occurs when the potential has an 
attractor as a response to a cosmological fluid, which ren- 
ders the final outcome (at fixed fiuid density) indepen- 
dent of the initial conditions. Many potentials support 
tracking behaviour, which takes place provided the initial 
dark energy density is not too low (otherwise the scalar 
field cannot dominate by the present). A tracking field 
slows down near the present as it starts to feel the friction 
induced by its own dominance of the energy density, and 
hence has equation of state reducing with time {w' < 0). 
Note that an asymptotic solution, in which radiation and 
matter are negligible and dark energy is the only compo- 
nent, may be an attractor but is not a tracker according 
to the above definition; tracking behaviour is induced by 
another fiuid component. 

A creeping field is one which sits static at low energy 
density until the density of other materials drops low 
enough for it to become dynamical and start to move at 
the present epoch. Creeping dark energy typically does 
not make predictions independent of initial conditions, 
but can be characterized by the equation of state increas- 
ing away from — 1 at the present epoch as the field starts 
to move [w' > 0). 

Potentials with tracking regimes also feature creeping 
behaviour for low enough initial density, which can lead 
to different predictions from the same potential. Creep- 
ers can be seen mainly as models which have passed the 
attractor and are frozen until they reach it, although 
creeping behaviour can also be found for many poten- 
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tials which do not support observable tracking.^ 

In order to regard the tachyon as a model of dynamical 
dark energy, it is useful to parametrize the barotropic 
index as a function of the scale factor: 



W = Wo+ Wa{a - 1), 



(13) 



so that w = today. Its value can be found by solving 
the Fricdmann equation today for Tq = T'(0), getting 
from Eq. (0 



= iUf - 1 = - 



1 - fiiQ 



From Eq. H10|l . one has that 

-N I 
Wa ^ e w , 

w' = 2wT' ( 3xT' + ^ 



(14) 

(15) 
(16) 



Then one has 



< when either T' < and 3xT' 



U^t/U < (for inverse power-law potentials, U ^ T ^ , 
the latter condition is 3xTT' - /3 < 0, true when T > 0), 
or T' > and 3xT' + U.t/U > 0. 

The transition between a tracker and a creeper can be 
defined by imposing an initial condition for T and T' so 



that 



0. Genuinely non-tracking models are prob- 



lematic as they do not make definite predictions, in the 
sense that there is strong dependence on initial condi- 
tions. This makes them rather unattractive, although 
by no means ruling them out. Indeed, it makes some of 
them more viable than trackers in relation to experimen- 
tal bounds, as we shall see. 



C. Choice for the tachyon potential 

As for the tachyon self-interaction, there are a number 
of models which one can consider, some being motivated 
by nonperturbative string theory and others purely by 
phenomenology. We review the classification of Ref. [3 , 
to which we refer for further details and references. For 
each case, past results are summarized. In the following, 
Ur is a normalization constant. 



I. U ^ UcT-P, /3 < 0, e cx T^-^ _^ qo as T 
This model is affected by instabilities p3 |. 



0. 



2.U^ UcT-f\ < /3 < 2, e cx T^^^ o as T ^ +oo 
[asymptotically de Sitter (dS)]. In order to get vi- 
able cosmologies, Uc does not need to be fine-tuned 



^ Whenever there is an attractor (either dust or de Sitter) in the 
phase-space plane, the asymptotic solution of course does not 
depend on the initial conditions. However, if these are 'too 
far away' from the attractor, it may not be possible to get 
close enough to a tracker/attractor before matter has reached 
its present density. 



since it is not affected by the super-Planckian prob- 
lem of the inverse quadratic potential. In general 
there is a stable late-time attractor . In Ref. 
the case (3—1 was considered as a numerical ex- 
ample. 

3. [/ = UcT~'^, e = constant. This is the potential 
associated to the exact power-law solution a — 
0,123, 13 ■ One has to fine-tune Uc in order to get 
sufficient acceleration today. A phase-space anal- 
ysis in Refs. ^, 25] confirms that it is difficult for 
this model to explain dark energy, since the only 
late-time attractor in the presence of matter or ra- 
diation has qt/x 0. 

4. [/ = UcT-/^, /3 > 2, e cx T^-^ ^ as T ^ +oo. 
This potential has not been studied numerically 
previously, but Ref. p3] showed that it has a dust 
attractor. The authors of Ref. ^ also argued that 
it behaves as in model 7 below. 

5. [/ = Ucexp{l/nT), ^ > 0, e cx exp{-l/ fiT)/T^ 

as T ^ +00. This gives an asymptotic dS solu- 
tion with effective cosmological constant given by 
Uc- This potential has been considered also for 
Klein-Gordon quintessence |0, l26||. and it should 
have properties similar to model 2 [9j. 

6. U = C/cexp(/i2r2), ^ > 0, e cx exp{- fi^T^) 
as T 0. This potential arises in KKLT setups 

for massive scalar modes on the Z?-brane j^l- 
The field oscillates around T — with effective 
cosmological constant given by Uc, and can give 
viable scenarios 9]. 

7. U = f/cexp(— /xT), /X > 0, e cx exp(/xT) ^ cx) as 
T +O0. This potential arises in D-D systems 
of coincident branes with a real tachyonic mode, 
and was studied in Ref. ■ The authors of Ref. 
found it to have a stable dust attractor after a pe- 
riod of acceleration. They suggested the possibil- 
ity that we are living during this transient regime. 
Note that this potential is a large-field approxima- 
tion of C/ C/c/cosh(^T) [23. 

8. [/ = Ucexp{-fi^T^), n> 0, e (X eyip{fi'^T^) ^ 
cx3 as T ^ +00. In Ref. @ it was argued that its 
predictions are similar to model 7. 

These potentials offer a rather comprehensive range of 
generic behaviours: one reaches a cosmological constant 
regime at infinite time [U cx exp(T~^)], one at finite time 
{U cx expT^), while in other situations the potential 
asymptotically vanishes. 

The structure of analytic asymptotic solutions of the 
equations of motion for inverse power-law and exponen- 
tial potentials is presented in the Appendix. 
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FIG. 1: Present value of the barotropic index w for tachyon 
tracking solutions (squares) and for creeping solutions for Ti — 
1 (circles) with inverse power-law potentials U oc T"'^. For 
comparison the result for tracking canonical quintessence is 
also shown (triangles). 



III. NUMERICAL ANALYSIS 

The equations of motion are integrated forward in 
time from before matter-radiation equality (iVj « —10).'* 
Having fixed the initial values Tj and T(, our numerical 
code adjusts the normalization of the potential to yield 
the solutions which satisfy the boundary conditions to- 
day (xq = 1, matter and radiation energy densities equal 
to flmfi and r^r.Oj respectively). In cases where a tracker 
behaviour is present, one would expect the present values 
of these to be independent from the initial conditions to 
a good approximation. 

We choose T( = as initial condition; an arbitrary 
non- vanishing value could typically only have arisen if at 
very early times the tachyon had an unacceptably high 
velocity. The initial phase of approach to a tracker is 
artificial, as the real cosmological initial conditions were 
presumably laid down at an earlier epoch than those of 
our code. Anyway, it will be sufficient to capture the 
main features of these models, although one might devise 
some physical situations which do predict U ^ g and 
< T' <C 1 as initial conditions. 

In the following we inspect case by case the models 
listed in Sec. Ill CI Since case 1 is unstable, we start with 
inverse power-law potentials. 



Although the problem is defined by boundary conditions at the 
present time, we do not attempt to integrate backwards from 
them (for instance by choosing Tq freely while Tg is constrained 
by the Friedmann equation at = 0). This is because in many 
cases we are dealing with situations having strong attractors, 
which become repellers in backwards integration leading to rapid 
growth of numerical instabilities. Instead, we integrate forwards 
using a shooting method to adjust the initial conditions to obtain 
the right present properties. 



A. U ^ UcT-'^ 

These potentials lead asymptotically to dS for /3 < 2 
and to a dust attractor for /? > 2. The first case is 
the most viable for obvious reasons, and we have verified 
that for Ti < 0.1 and T( = the numerical solutions with 
/3 < 2 lie on the asymptotically dS attractor, while for 
(3 > 2 and the same initial conditions it is not possible to 
achieve a cosmology compatible with observations (i.e., 
xq = 1). For Ti > 0.1, one goes further and further away 
from the attractor and enters a creeping regime, in the 
sense that the evolution depends on the choice for 7]. In 
Fig.nthe value of wq is shown for the tracker solution and 
for a particular creeping solution with Ti — I. In general, 
creepers mimic a cosmological constant more closely than 
trackers, and, in particular, only a creeping solution is 
available for (3 > 2. For comparison we also show the 
equation of state from tracking canonical quintessence 
models,^ which for the same /? is always more negative 
than that of the tracking tachyons. 

These results can be understood by recalling that in 
the slow-roll approximation there is a map pin g between 
the DBI scalar and the Klein-Gordon one [mTisII . The 

tachyon energy density is gr = U{T)/\/l - U{T ) + 

U{T)f^/2, and after a field redefinition cj) = J ./U{T)dT 
one obtains a canonical action for a scalar (j) with poten- 
tial W{(l)) = U[T{(t))]. For U{T) - T"/^, one has 

W{<i^)^rP , /3 = ^- (17) 

Hence the tachyon tracker for < /3 < 2 is dual to the 
tracking quintessence solution /3 > 0. From Fig. ^ one 
can see that any dual pair (/3, (3) corresponds to almost 
the same index — see in particular the pairs (2/3, 1) 
and (1,2). 

The attractor evolution of the fluid components is 
shown in Fig. [5] together with w{N), for the particular 
case /3 = 1/3. The behaviours of these quantities for 
other values of /3 are all similar. Note that the barotropic 
index is w ~ —1 until matter domination, increases up to 
w « —0.8 at redshift z ~ 10^, and then decreases towards 
— 1. This is consistent with the 'freezing' behaviour and 
its limit of applicability as classified by Ref. [20l |. 

For the creeping solutions (not shown here) the index 
w is very close to —1 up to very late times, when it starts 
deviating to a softer equation of state. 

In this and all other examples, the onset of dark energy 
lies in the interval — 1 < iV < 0, that is, for redshift z < 
1.7. Such late domination by dark energy is essential to 
prevent an excessive suppression of structure formation 
growth, see e.g. Ref. 0- 



^ Canonical quintessence follows a Klein— Gordon equation (j> + 
3Hip + V^,f, = 0. The tracking regime is reached from small initial 
values of the scalar field ij> (</>; <C 1 in Planck mass units). 
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FIG. 2: Evolution of the density parameters and equation 
of state as a function of A'^ for the inverse power-law tachyon 
model with jS = 1/3. Solid line: g-m/x; dashed line: Qr/x\ 
dotted line: qt /x. The dot-dashed line is the barotropic index 
w. 
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FIG. 3: Numerical points in the wq-Wq plane for tracking 
solutions (squares) and creeping solutions for Ti — 1 (circles) 
with inverse power- law potentials U oc T~^. The value of (3 
is shown at each point, while the dotted lines are the la and 
2a likelihood bounds llal . The unlabelled points on the left 
and right of /3 = 1 in the tracking curve are f3 = 0.95 and 
l3 = 1.05, respectively. 

In Fig.Othe points predicted by the same models in the 
wq-w'q plane are shown together with the la and 2a like- 
lihood contour bounds of Ref . . These are based on 
the first-year SNLS data set [H and SDSS [13, WMAP3 
[T^ . and 2dF [^113 experiments. Note that there is a 
minimum value for w'q around P « 0.99; we have checked 
that its closeness to 1 is an accident of the particular 
f^m,o used. 

The conclusion is that there is tracking behaviour for 
/3 < 2, and if Ti is not too large this operates and gives 
a well-defined prediction for wq, which is however in dis- 
agreement with observations unless /3 is quite small. The 



predicted wq moves closer to —1 as the potential flat- 
tens, but only once it gets down to about /3 < 0.8 do 
we start to see compatibility with the 2cr-level bounds. 
This is even flatter than the equivalent constraint on 
quintessence, f3 ^ 4/3. Hence these models are not very 
pleasant from either a theoretical or experimental point 
of view. One has to impose /3 ^ 1 or fine-tune the initial 
conditions. Relaxing either assumption results in incom- 
patibility with observations. The creeping solution with 
Ti — 1 (upper curve) gives an example of observation- 
ally viable model for /3 < 3, but at the cost of fine-tuned 
initial conditions. 

One can also make T( big enough that the initial w is 
close to zero rather than —1. Then it dips to —1 and 
increases to join the tracker, eventually giving exactly 
the same evolution as if one had chosen 7]' = at the 
start. This is exactly the same sequence of epochs that 
are seen in quintessence models '33] ; if the initial velocity 
is set high, then the trajectory 'overshoots' the tracker, 
comes to rest at lower potential energy for a while, and 
eventually rejoins the tracker from below. 

We can compare the numerical output for Tq with the 
asymptotic solutions derived in the Appendix, specifi- 
cally Eq. PTejl . With C = 0.045, the ratio |ro"""> - 
ytheor|/y„um jg ^Iways Icss than 2%. As explained in the 
Appendix, the asymptotic analytic solution relies on a 
polynomial rather than monomial potential, and in gen- 
eral they will predict a different variation oiw. Moreover, 
the values of p from Eq. (|A15|) are rather small, corre- 
sponding to —2.5 ^ N^, ^ —0.5, where N^, is the big bang 
event; i.e. the asymptotic solution which matches the 
present evolution has a big bang in the very recent past, 
and hence can only have become a reasonable approxi- 
mation very recently. The solution properly approaches 
the attractor only for large positive N, when the con- 
tributions of radiation and dust matter are negligible. 
Nevertheless it appears to give a good estimate of wq, 
though not Wq. 



B. [/ = I7cexp(l//ir) 

This model has a non-zero vacuum energy Uc as 
T — > oo, which will be achieved asymptotically. It ex- 
hibits some features similar to the inverse power-law po- 
tential with (3 < 2, but does not appear to exhibit track- 
ing behaviour. We have fixed fJ, = 1 and checked that 
this does not result in a loss of generality. For Tj > 2 
the solution is very close to de Sitter, wq ~ — 1 and 
< Wq < 10~^, essentially amounting to a creeping solu- 
tion [Tq « Ti) sitting on the asymptotic flat part of the 
potential. 

For smaller Ti < 1 the cosmological values today de- 
pend on the initial conditions of the model. For Ti = 1 we 
find Wo ~ -0.95, w'o « 0.04 and for T; = 0.5, wo « -0.77, 
Wq w —0.05. The solution exits the 2a contour at about 
Ti = 0.4 {wo « -0.70, w'q w -0.14). One might ar- 
gue that the initial condition T( = is inappropriate for 
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this model, since the potential is very steep at small field 
values. However, we have verified that the slow-roll ve- 
locity is small enough to be negligible anyway at least for 
Tj > 0.2 (which is beyond the observable region). There 
is no evidence of tracker- type behaviour for Ti > 0.2. 

This model does not need fine-tuning of the initial con- 
ditions, as the field can lie anywhere on the flat part of 
the potential. However one has to fix the value of Uc 
by hand in order to reproduce the observed dark energy 
density, which of course does not resolve the cosmological 
constant problem. 

As in the inverse-power case, we can compare the 
numerical output of wq with the asymptotic solution 
Eq. (IX35ll . where </>o = i/ynum y^ ^ien 0o - 1 
and the approximation giving Eq. (jA35p breaks down, 
ko"™ - ~ 0-1' while at large Tq one has 

j^nu m thcorj^j^^numj < ;^q-3 q^. better. Wc UOtC that 

Eq. ljA36l) is negative definite while wg"™' > 0, and there 
is disagreement between the two. 

C. U ^Ucexp{^l^T^) 
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FIG. 4: Evolution as a function of A'^ of the model U = 
Uc/ cosh(T), in this case extending beyond the present. Solid 
line: Qm/x; dashed line: Qr/x; dotted line: qt/x. The dot- 
dashed and dot-dot-dashed lines are w and w' , respectively. 
In this example, the initial condition is Ti — 0.5, for which 
wo ^ -0.98 and w'o ^ 0.03. 



This model also has a non-zero vacuum energy Uc, this 
time at T = 0. Fixing /j, = 1, for Tj < 0.1 one has a de 
Sitter behaviour {wq w —1, < Wq < 10"'^), while for 
increasing T; the barotropic index goes away from —1 (for 
instance, Ti = 1 gives wq ~ —0.81, w'q « 0.13). 

Checking the numerical output of wq with Eq. ljA37|l . 
where (jjo = (Tq""™)^ < 1, one can see th at jwg"™ - 
^theor|/|^„um| < ;^q-4 Again, Eq. |J38|) is negative 
definite while wg"""' > 0. 

Note that lower values of fi lead even more closely 
to a cosmological constant behaviour, well inside the la 
bound. This may be relevant when trying to construct a 
model which fits also for infiation, as one needs fi^ ~ 10^^ 
to get the correct level of anisotropics . 



dust regime are not very predictive. They are however 
capable of explaining the observed dark energy proper- 
ties. 

E. U ^Ucexp{-ii^T^) 

This model has a dust attractor whose qualitative fea- 
tures match the previous case. It is not difficult to find 
suitable initial conditions mimicking a cosmological con- 
stant today. 

IV. DISCUSSION 



D. U ^Uc exp(-^iT) and U = Uc/ cos]i{fiT) 

When considering the pure exponential potential U = 
Uc exp{—fiT), choosing a different starting value of T has 
no impact on the evolution because a rescaling T — > T+ 
const, simply renormalizes Uc, which the program then 
adjusts to give the same present status (corresponding to 
Wo « —0.93, Wq ~ 0.10). Therefore we consider directly 
the hyperbolic cosine potential of Ref . , for which this 
degeneracy is removed. For < Tj < 0.3, the solution has 
Wq ^ —0.99 and < Wq ^ 10~^, while for larger values of 
Ti the present-day barotropic index becomes larger. We 
note that the accelerating phase of these solutions is only 
a transitory epoch before reaching the dust attractor in 
the future, as shown in Fig. 0] where the equations have 
been integrated up to positive values of N. 

Since one must tune, although not too severely, the 
initial condition (or, equivalently, Uc) so as to get viable 
acceleration today, this and other models leading to a 



None of the models we have discussed are very satisfac- 
tory. Those which carry reasonable theoretical motiva- 
tion all end up with a high degree of fine-tuning. Either 
the potential normalization Uc has to be set to match 
the observed dark energy density, or the initial condi- 
tions tuned to the creeping regime, meaning that the 
present density was already set in place during the early 
Universe. In either case, this tuning amounts to a re- 
statement of the cosmological constant problem, rather 
than a resolution. The models which avoid fine-tuning 
of initial conditions (while still subject to tuning of the 
normalization), such as the inverse power-law, are con- 
strained by observations into parameter regimes with no 
theoretical motivation. 

The string tachyon with DBI action seems only weakly 
competitive as a dark energy candidate, and another for- 
mulation of the effective theory might have more success- 
ful applications. The tachyonic effective action as the 
lowest-order level truncation of cubic string field theory 
has been studied only very recently and its cosmological 
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impact has yet to be fully assessed [s^, 'ss', '36', . An al- 
ternative would be to abandon the pure four-dimensional 
picture (ideally corresponding to a low-energy single or 
coincident brane configuration in a higher-dimensional 
spacetime) and consider more general brane-anti-brane 
setups where open string modes naturally live. However, 
any modification to Einstein gravity would be severely 
constrained after nucleosynthesis. 

Finally, it was shown that the addition of electric and 
magnetic fields in the non-BPS brane world volume can 
slow down the evolution of the infiationary tachyon and 
relax the fine tuning of the parameters [33, |33 ■ It might 
be worth checking whether the insertion of non-trivial 
fluxes plays some role in the late universe. 

Still, the era of high-precision cosmology opened up by 
microwave background and large-scale structure observa- 
tions is allowing us to constrain infiationary and dark en- 
ergy models in a more and more stringent way, selecting 
some of them from a plethora of possibilities. Here we 
have given an example of the first stage of such a proce- 
dure. The second one will be to develop and refine those 
models that seem particularly promising, by embedding 
them in a comprehensive and consistent picture of the 
cosmological history and its particle theory content. 
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APPENDIX A: ASYMPTOTIC SOLUTIONS 

In this Appendix we derive some analytical asymptotic 
solutions describing the evolution once the tachyon has 
become completely dominant. 

1. Asymptotic solutions for power-law potentials 

Inverse power-law potentials, introduced in the context 
of dark energy in Ref . |40l | , have no general support from 
string theory. However, they can be viewed as large-field 
approximations (|T| — > -|-oo) of the exact solution 

a = exp[p{t/toT -p], (Al) 
where we have normalized a so that oq = 1. Wc have 



2(2-/3) 



(A4) 



for < n < 2, or 

U ^ — rp-i3n-i)/{n-2) 

3-7 



(A5) 
(A6) 



for either n < or n > 2. The solutions with the poten- 
tial Eq. ljA3|l are real, nontrivial and expanding if, and 
only if, 

iV, < and < /? < 2 (0 < n < 1) , (A7) 



iV, > and /3 < {1< n < 2) , (A8) 

where iV* = —p. In the first case, the big bang event is 
at iV* and N^, < N < +oo; in the second case, — oo < 
N < N^,. The solutions with the potential Eq. ljA5|l are 
well behaved when iV* > and either 2 < 7 < 3 (n < 0) 
or 7 > 3 (n > 2). 

The case given by Eq. (jA7|l is the most interesting 
since it corresponds to the tracking regime. Using the 
definition Eq. l|A4p and neglecting matter and radiation 
contributions, one has the exact solution 



x 



w = l = (wo + l)-r^ 77 

3(2-/3) , 2 
w = (w + 1) , 

y/2 ^ + 1 



(A9) 
1, (AlO) 
(All) 



T = 



(wo + 1) 



2(1-/5)7(2-/3) 



N,-N 



(A12) 



Sy/wo + 1 



1/(2-/3) 



C, (A13) 



where we have assumed T' > (T > 0), C is an integra- 
tion constant, and 



wo + l = - - 







1 



3(2 - /3) N, 



(A14) 



One can find the value ofp for the attractor from today's 
value of the barotropic index. Inverting Eq. (jA14|) . 
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Also, 



To 



3(2-/3)(w;o + l)' 
f3 



3V'iA'o + 1 



Tq = Vwo + 1 , 
3(2-/3) 



Wn 



13 



+ C, 



(wo + If. 



(A15) 

(A16) 
(A17) 
(A18) 
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The classical stability of these solutions was studied in 
Refs. 0, At late times (i.e., large T and N), one 
has a de Sitter regime for < /3 < 2, e 0, or an 
asymptotically dust solution (7 > 2). These solutions 
are actually attractors. As regards the latter, we note 
that 



T2 ^ + 1) 



1, 



(A19) 



where we have taken the Carroll limit (the DBI ac- 
tion then becomes singular) corresponding to tachyon 
condensation into dust. This expression prescribes a reg- 
ularization for all the above formulae, and suggests the 
following possibility, considered also in Ref. . The only 
way to balance the increasingly small denominator is to 
impose that wo + 1 ~ 0; this might mean that if 7 > 2, 
the Universe tends to become dust dominated but only 
after passing through an accelerating phase. Since the 
origin of time is arbitrary, such a phase is not positioned 
unequivocally and it will be determined also by the nor- 
malization constant Uc of the potential. 



2. Asymptotic solutions for exponential potentials 

In general, a solution for the Friedmann equation in 
terms of T can be found by noting that, when w ~ const, 
the T dependence of the Hubble parameter must be the 
same as of the tachyon potential, the square root in the 
denominator of Eq. (j^ (with g — 0) being dimensionless: 
x{T) cx U{T). All the exponential potentials we consider 
can be suitably parametrized so that 



X = Aexp[(^T)"] , 



(A20) 



where A = exp[— (/iTo)"] and we have neglected all mat- 
ter/radiation contributions. Differentiating this equa- 
tion with respect to N and using the continuity equation 
{w' = 0), one has 



which can be integrated from N to today: 



(A21) 



N 



— / drri-"exp[(A*T)"]. (A22) 



Defining the variable = (fiT)", the integral becomes 



iV 



3A 

3^(_l)2/a 



/j2(l-a)/ag0 



(A23) 



r --1,-00 

a 



Fia, 



(A24) 



where T is the incomplete gamma function. The cases of 
interest are 



fe*^"*^" - 1) 



(A25) 



F{-i,ct^) = ^^[r(-3,-0o)-r(-3,-0)] 



F(2, 



3e- 



The first equation implies that 



(A26) 

[r(o,-0)-r(o,-0o)]. (A27) 



T{N) = To + j^ In (^1 - 



(A28) 



In order to find T as a function of N in the other cases, 
Eqs. l|X26|) and l|X27)l can be expanded around large or 
small d: for a — ~1, 



«l: F(-1,<A) ^ ^(^-^).(A29) 



(A30) 



while when a — 2, 



n.,„._(_-_).,A32) 

One can then numerically invert the above transcenden- 
tal functions to get (l>{N) and 



w = -1- (t)'{N)/3, 
w' = -</)"(7V)/3. 



(A33) 
(A34) 



When (p ^ 1 (which is the regime typical of the dS- 
attractor solutions), one can show that 



33 

for a = — 1, while 



o ^J. e (Pq , 



2/i 



wo ^ 1 Y ) e^"(j)o - 1 , 



33 



(A35) 
(A36) 

(A37) 
(A38) 



for a = 2. Note that w' is always very small in these 
models. 
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3. Comparison with numerical solutions 

Having built the asymptotic analytic solutions, one 
would like to check whether the numerical behaviour re- 
ally approaches such solutions at late times. In partic- 
ular, we want to compare the numerical points in the 
Wq-w'q plane in parameter space, found via 



Wo 



Wn 



2woU 



TO 



(A39) 
(A40) 



with the corresponding semi-analytic expressions given 
by Eqs. (|IT6)l-ljIT8)l fo r inverse power-law potentials 
and by Eqs. ljA35p - (|A38(l for exponential potentials. The 
matter content today still amounts to 25% of the total 
energy density and contributes to the cosmological evo- 
lution, so we expect a deviation of the numerical results 



relative to the asymptotic, pure tachyonic solution. 

There is another source of discrepancy one should take 
into account, namely the approximations implicit in the 
solutions presented in this Appendix. In the inverse 
power-law example, while Eqs. ljA9|l - l|A13|l are valid for a 
polynomial potential as in Eq. (|A2(I , the numerical model 
is actually Eq. (|A3|) . and in general the two will give a 
different running of the barotropic index w' , Eq. ljA40p . 
In Sec Him we find that Wq""™ and Wq do disagree, 
but still there is remarkable agreement between w™™ 
and wji^™'' [the latter comparison is in fact done between 
To"""" and r(5''<="(u;S"™)]. 

In the case of exponential potentials, the asymptotic 
models have w « const, by construction [x{N) cx U{N)] 
and the equations of motion are much simplified. Again, 
^num ^ ^thoor good approximation (see Secs. lIII Bl and 
nil C|) . while a comparison of the running would require 
a more refined treatment. 
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